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1. Introduction 

A dimensionally reduced effective theory is, by now, a well-established and powerful tool 
for studying high-temperature QCD both perturbatively and non-perturbatively. At high enough 
temperatures the QCD coupling becomes small, and perturbative methods can be safely applied to 
hard (~ T) modes. The fact that high-temperature QCD can be reduced to 3d effective theory can 
be understood by considering an Euclidean finite-r propagator 

1 

p 2 + CO 2 + m 2 ' 

where 03% = 2nnT for bosons and col = (2n + \)%T for fermions. Thus, only static n = bosonic 
modes remain "light" at high temperatures, and the heavy n ^ bosonic and all fermionic modes 
can be integrated perturbatively [jl], 0, ^, The result is an effective 3d theory of full QCD, 
electrostatic QCD (EQCD) JIT|]. An essential feature of the effective theory is that it fully includes 
the perturbatively problematic soft (~ gT) and ultrasoft (~ g 2 T) scales of the original theory; only 
perturbatively well-controlled hard scales (~ T) are integrated over. 

EQCD offers a good starting point both for perturbative calculations []5j,|6|, [7|] and non-perturbative 
lattice simulations. In the latter case, EQCD offers an interesting alternative to standard high- 
temperature lattice simulations: above all, the effective theory is purely bosonic and only 3- 
dimensional, making it much cheaper to simulate. The theory is superrenormalizable, which ren- 
ders the continuum limit particularly transparent; it also enables simulations at arbitrarily large 
temperatures. On the other hand, the effective theory cannot be used to study the QCD phase tran- 
sition: at too low T QCD becomes strongly coupled and the perturbative derivation of EQCD fails. 
Nevertheless, the theory has been found to work well down to temperatures T ~ 1 .5 — 3 T c , depend- 
ing on observable used. Lattice simulations of EQCD have been used to calculate QCD pressure 
at high T [j|], spatial string tension []9|], and spatial screening lengths JTOj]. Here we measure the 
quark number susceptibility, and compare the result to both perturbative and 4d lattice results. 

2. Susceptibility in electrostatic QCD 

EQCD is defined by the action 

S E = J d\^E 

&l = ^Tr[i$ + Tr[D,-,A ] 2 + m^Tr^] + i^Tr\A 3 ] + ^(Tr^]) 2 , (2.1) 

where Fjj = djAj — djAj + igi[Aj,Aj] and D, = dj + igjAj. Fjj, A,- and Ao are traceless 3x3 Her- 
mitean matrices (Ao = A^T a , etc). Coupling and mass parameters g$, m^, 73 and A3 are defined by 
the physical 4d temperature, renormalization scale Aj^-, chemical potential \i and the number of 
massless fermions. It is convenient to use the dimensionless ratios 



Quark number susceptibility at high temperature 



Ari Hietanen 



which determine the physical properties of EQCD. The pi -dependence of the parameters is, at 1- 
loop level, 



where pi = pi/(nT) and the pi =0 expressions can be found in ref. [11]. The two loop corrections 



have been calculated in ref. [10], but the effects remain in practice negligible. 

The quantity we are interested in is the quark number susceptibility, which we calculate over 
one flavor u only, can be defined in EQCD as 

%3 = v^ ln ^ = v~^ ln .f mkAo exp (-5e) (2 - 4) 

Substituting Se from ( |2~l| ) we arrive at equation 

^3 = -2^^=o(T^) + y§((T^) 2 ) 
9N 2 f 

+V (2N C +N)^ 2 rt=° ( ((TrA ° )2) " ( ™« )2 ) (2 " 5) 

Thus, the quark number susceptibility is obtained by measuring the condensates (TrA^), ((TrA 2 ,) 2 } 
and ((TrA 3 ,) 2 ) on the lattice. Due to the superrenormalizable nature of the theory, measurements 
can be rigorously converted to MS scheme in the lattice continuum limit; because MS was used in 
in the perturbative matching to 4d QCD, this also allows us to compare to 4d results. For example, 
the continuum limit for the ([7YA 3 ,] 2 ) is 

V<[TrA 3 /g 3 ] 2 >^ = hm M |y([TrA^]2]) a __A_ [lni3 + o.08848010]|, (2.6) 

where j8 = =^ is the lattice coupling constant. Continuum limits for (TrAfy and ((TrA 2 ,) 2 ) are 
given in [O]. The relation between %}, an d the ttue 4d susceptibility is 

* = f^ + « A "' <2 - 7) 

where Ap = pqcd — P3d is the perturbative 3d— >4d matching coefficient, and can be found in [§]. 

3. Lattice measurements 

Lattice simulations were carried out for Nf = 2. We used large lattice sizes V = 140 3 ,200 3 and 
values of j8 =32, 40, 67, 80, 120. The volumes are large enough so that the remaining (exponential) 
finite volume effects are safely below our statistical accuracy; for a detailed analysis in a related 



theory see [12]. 

For each value of /3 we did simulations with values of y = 6.62, 5.31, 4.00, 3,09, 2.02, 1.18, 
0.71, and 0.45. The simulations have been carried out for pi = 0, and to obtain the susceptibility, 
one needs measurements of (TrA 2 ,) and ((TrAg) 2 }. 



3 



Quark number susceptibility at high temperature 



Ari Hietanen 



y=6.62 



y=5.31 



y=3.99 




H A 

CO 



-0.76 



^n^t £ -0.8 — 




-0.82 — 



0.01 0.02 0.03 
1/P 



0.01 0.02 0.03 
1/P 



0.01 0.02 0.03 
1/P 



% =40.5/3 
5C 2 =6.1/2 



X =13.3/3 
%= 1.53/2 



X =11.9/3 
5C 2 = 1.82/2 



Figure 1: Difference of continuum extrapolations with and without the log-term. The solid line is a fit 
Cl + J + TP and the dashed line c i + f + J l °g(P) - 



C3 



The standard procedure of doing continuum extrapolation is to fit a polynomial to the diver- 
gence subtracted lattice data. However, in EQCD the divergent lattice contributions contain also 
terms of the form j3 log(/3) and such terms could also arise in terms ~ a ~ 1//3. Therefore, we do 
continuum extrapolation in two ways; we use a second order polynomial 



and a fitting function of the form 



ci + J + |log(/3) + ||. (3.2) 

The Fig. [T] demonstrates the differences between different continuum extrapolations for (TrAg). 
The extrapolations done using (3.2) have excellent # 2 /dof, and we shall use this form henceforth. 



The full continuum extrapolations are shown in Fig. g. We note that the detailed form of the 
fitting function is significant for our results; thus, knowledge of the true l/j8 coefficient is highly 
desirable. There is an ongoing calculation of this term using stochastic perturbation theory [jl3|], 
which will hopefully confirm our result. 

The contributions of ((TMq) 2 } to susceptibility are numerically much smaller than that of 
(TMq}. The measurements of ((TrA ( 3 } ) 2 ) are not accurate enough to distinguish between extrapola- 
tions (3. 1) and (|3T2|), and the choice is not significant for the final results. Hence we use the second 



order function in Fig. |2| 

The continuum extrapolated results are shown in Fig. ||, as a function of the parameter y(T). 



The result agrees very well with the perturbative susceptibility derived in ref. []14J. The perturbative 
result has the form: 

Zpert = a i y i/1 + a 2 y + a 3 y 1/2 + a 4 (3.3) 
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Figure 3: Le/L - the susceptibility as a function of y and T/A-^. The solid line is the perturbative prediction. 
Right: the difference between the perturbation theory and lattice results. Points denoted by squares are 
obtained from logarithmic continuum extrapolation and diamonds from polynomial extrapolation; the former 
one has expected behavior as y — > °°. 



where a\ are functions of x and z. On the right panel we show the difference between the per- 
turbative result and the lattice measurement. Because the lattice measurement (in the continuum 
limit) fully includes perturbative contributions, we expect the difference to behave as 0(y^ 1 ^ 2 ). 



This is the case when continuum extrapolation is performed using Eq. ( |3.2[ ); the extrapolation ( p.ip 
diverges from perturbative result as y — > °°, which is inconsistent. This also motivates the use of 
(^2) as continuum extrapolation. 

In Fig. [| we show the susceptibility in 4d units, transformed using Eq. (2.7), and compare the 
result to 4d Nf = 2 staggered fermion lattice simulations [If], 17]. Here we use the ratio T c /Ay^ = 
0.49 [ 1 8 ] . We can observe that our results agree very well with 4d simulations and with perturbation 
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Figure 4: The quark number susceptibility in 4d units. Our results are shown with diamonds, only 3 points 
fit in the temperature range shown. The results agree very well with 4d lattice simulation results JjjL Il7j], 



theory. Indeed, the difference between EQCD measurements and perturbative results are all but 
invisible on the scale of the plot, possibly excluding the lowest temperature data. 



4. Conclusions and Outlook 



We have measured the quark number susceptibility in dimensionally reduced effective theory 
of high-r QCD. The results show a good agreement with the perturbation theory and 4d lattice 
simulations, down to temperatures ~ 3T C . Currently we are expanding our simulations for finite 
chemical potential. In order to avoid the sign problem at finite jJ. we do the simulations using 
imaginary chemical potential and then continue analytically to real values. 
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